Spiking Optical Patterns and Synchronization 
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We analyze the time resolved spike statistics of a solitary and two mutually interacting chaotic 
semiconductor lasers whose chaos is characterized by apparently random, short intensity spikes. 
Repulsion between two successive spikes is observed, resulting in a refractory period which is largest 
at laser threshold. For time intervals between spikes greater than the refractory period, the distribu- 
tion of the intervals follows a Poisson distribution. The spiking pattern is highly periodic over time 
windows corresponding to the optical length of the external cavity, with a slow change of the spiking 
pattern as time increases. When zero-lag synchronization between the two lasers is established, the 
statistics of the nearly perfectly matched spikes are not altered. The similarity of these features 
to those found in complex interacting neural networks, suggests the use of laser systems as simpler 
physical models for neural networks. 



PACS numbers: 05.45.Xt, 42.65.Sf, 42.55.Px 

Semiconductor lasers, subjected to optical feedback, 
display chaotic behavior The chaotic behavior con- 
sists of a very short and random spiking of the laser inten- 
sity with the time between spikes depending on how far 
above lasing threshold the laser is. Two chaotic lasers 
can be synchronized with each other and this has al- 
lowed them to be excellent candidates for novel broad- 
band [3, H, 0, [E 01 communication devices. Different 
configurations, such as delayed optoelectronic 0, [1| or 
coherent optical injection 0, [13, [HI have been used 
for synchronization of the two lasers. Using optical feed- 
back, configurations consistingof unidirectional 0, Q or 
mutual coupling @, |^, l9l.Tl0l. [ll| and variations of 
the strength of the self and coupling feedback have been 
shown to result in different synchronization states. The 
lasers can synchronize in a leader-laggard or anticipated 
mode, as well as in two different synchronization states; 
achronal or generalized synchronization [13) [13) [13] where 
the cross correlation is time shifted by the feedback de- 
lay time but neither laser acts as a preferred leader or 
laggard, or isochronal synchronization (zero-lag) where 
there is no time delay between the two lasers' chaotic 
signals |,ilil,[il,|i3. 

Zero-lag synchronization of lasers was recently ex- 
tended to a cluster consisting of three semiconductor 
lasers, mutually coupled along a line, in such a way that 
the central laser element acts as a relay of the dynam- 
ics between the outer elements [3) [3l- The zero-lag 
synchronized dynamics of remotely located chaotic signal 
sources has sparked an interest in such systems in part 
because they have features also seen in biological and 
neural transmission networks. Though the time scales 
for the two phenomena are vastly different; lasers spik- 
ing on 100 ps time scale while neurons spike on ms time 
scales, much of the dynamics and spiking statistics ap- 
pear to have common behavior. Here we report on the 
spiking optical pattern of solitary and two mutually cou- 
pled chaotic lasers, observed on a time scale which re- 
solves the individual spikes, in both their synchronized 



Laser A 


-A ►— 




Laser B 




— < ►- 



partially transmitting mirror 



FIG. 1: Schematic experimental setup. Two lasers are mutu- 
ally coupled via a partially transmitting mirror placed in the 
middle of the coupling optical path between the lasers. 

and unsynchronized states and determine the statistical 
behavior of the spiking. This further allows us to es- 
tablish an analogy to the spiking behavior of single and 
interacting neurons. 

Our experimental setup is shown schematically in Fig- 
ure 1, where two semiconductor lasers are coupled via a 
partially transmitting mirror placed in the middle of the 
coupling optical path between the lasers. In the actual 
experiment the self feedback and the mutual feedback 
paths were spatially separated through the use of beam 
splitters The time for light propagation from the 
mirror to one of the lasers is t/2 . We distinguish be- 
tween the following three limiting scenarios. In the case 
of a fully reflecting mirror, the lasers are uncoupled and 
each laser is subject only to a delayed self-feedback. In 
the second scenario where the mirror is fully transparent, 
each laser receives a delayed signal from the other, and 
this configuration is known as " face-to- face" . In the in- 
termediate scenario, the mirror is partially transmitting 
and each laser is driven by two delayed signals, one from 
self-coupling and one from mutual coupling. 

In the experiment we used Fabry-Perot semiconductor 
lasers emitting at 670 nm wavelength, selected to have 
nearly the same threshold current, emission wavelength, 
and output power. The temperature of each laser is stabi- 
lized to better than Q.QIK and the individual laser tem- 
peratures are tuned so that the two lasers have nearly 
identical output wavelengths. The self and mutual feed- 
back loop time is 23.55 ns. Two fast (50 GHz bandwidth) 
detectors biased via a 40 GHz bandwidth bias T measure 



the output intensity of each laser. The DC current into 
the bias T is used to measure the average DC power 
faUing on the detector while the AC currents are mea- 
sured simultaneously in two channels of a 12 GHz band- 
width, 40 GbGS/s oscilloscope (Tektronix TDS 6124C). 

For the case in which the mirror is fully reflecting the 
lasers are decoupled. Each laser becomes chaotic due to 
the self-feedback but their chaotic fluctuations are com- 
pletely independent of each other. A typical trace of one 
of the laser output intensity measurements is presented 
in the top panel of Figure 2a, where the ratio of the ac- 
tual laser current to the threshold current, p, is 1.03. The 
time dependent intensity of the laser consists of spikes of 
typical duration of ~ 120 ps. In the following panels we 
show the same laser intensity fluctuations recorded af- 
ter a time, t(= 23.55 ns) , 2t and 3t. Figure 2 clearly 
indicates that on a time scale of a few optical delayed 
self-feedback times the timing of the spikes repeats itself, 
where as time elapses the spikes gradually broaden and 
finally disappear as new spikes emerge in new positions, 
forming a new pattern. This behavior is physically easy 
to understand since the feedback photons circulate in the 
long external cavity with a periodicity r and the lasers re- 
ceives nearly the same feedback pattern with this period- 
icity. The photon lifetime in the cavity is finite, however, 
and thus the feedback waveform slowly changes. After 
many r periods (many photon round trips) the feedback 
waveform and the laser's chaotic fiuctuation pattern will 
have changed completely. From this explanation we can 
also see that the revival of the chaotic pattern with period 
T will be nearly independent of p since the chaos is caused 
and determined by the feedback photon train which has 
a similar form with period r. To further demonstrate the 
repeatability of the patterns after a delay t in Figure 2b 
we show the overlaid intensity trace at time t-\-T (green) 
and at time t + 2T (blue) . 

A quantitative analysis of the spike statistics requires 
the definition of a low-signal threshold so as to eliminate 
the small spikes in the measured laser intensity which 
are the result of noise in the measurements. Our thresh- 
old was chosen to be at twice the average detector noise 
level, and the timing of a spike above this threshold level 
was determined according to the time of its maximum 
intensity. For each measurement we accumulated 70, 000 
consecutive spikes. A histogram of the time interval be- 
tween consecutive spikes is presented in Figure 3a for 
various values of p. 

The probability distribution of the time intervals be- 
tween two spikes consists of two main features. The dis- 
tribution for relatively long time intervals between the 
spikes follows a random, Poisson distribution where for 
small p values the exponential decay rate increases lin- 
early from zero as a function of p, as shown in Figure 
3b4a. For very short time intervals the Poisson distribu- 
tion is altered so that immediately after a spike it is most 
improbable to record a second spike. The most prob- 
able time between consecutive spikes is defined, as for 
neural spike trains, as the refractory time [13, [2l[. Fig- 




FIG. 2: (a) A trace of 15 ns duration of the intensity of one 
laser followed by plots of the same laser intensity after a time 
r, 2r and 3r with r = 23.55 ns. The laser was operating with 
p — 1.03 and with a reflected power of a few % of the laser 
output intensity, (b) The intensity trace at time t + r (green) 
and at time t + 2t (blue) demonstrating the slowly decaying 
periodicity of the spiking pattern. 



ure 3c 4b indicates that the refractory time increases as 
the laser current approaches the threshold value, which 
makes physical sense, since at low pumping currents it 
takes longer to rebuild the laser gain after the previous 
pulse had depleted it. 

The chaotic dynamics are well described by the three 
coupled Lang-Kobayashi rate equations for the optical 
field amplitude E, the optical phase $ and the excited 
carriers, n, of the gain medium [2^. The equations also 
predict that the intensity chaos is due to a spiking be- 
havior of the lasers on a 100 ps time scale, with no laser 
emission between spikes, in full agreement with the obser- 
vations. Simulations From the simulations of the Lang- 
Kobayashi equations indicate we also determined that 
after a spike the laser ceases to emit because the number 
of carriers drops well below threshold, hence a successive 
spike is forbidden until the population is repumped again 
by the laser injection current. 

For the case of a fully transparent mirror, the two 
face-to-face lasers synchronize achronally, characterized 
by two dominant peaks in the cross correlation function 
at ±r. We find that the statistics of the intervals between 
spikes, the refractory period and the repetitive spiking 
pattern with optical delay time, r, are not altered. How- 
ever, when When the two face-to-face lasers are config- 
ured to have different p values, the lasers no longer syn- 
chronize, and from a sample of the signals transmitted 
between the lasers it is not obvious that one could not 
determine that the two lasers are operating with two dif- 
ferent p's and the values of the two p's. the The statistics 
of the spike intervals, however, is altered. reveals tis in- 
formation. Figure 4 5 depicts the distribution of times 
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FIG. 3: The probability for time intervals between spikes for 
various values of p. 




FIG. 4: (a) The exponential decay rate as a function of p ob- 
tained for time intervals greater than the most probable time. 
The decay rate is a linear function of p for low p values, as 
indicated by the linear fit (solid line), (b) The most probable 
time interval, the refractory time, as a function of p. 



between spikes when pi ~ 1.1 and p2 ~ 1.4 where as 
previously, the feedback strength for each laser is a few 
percent of it's output power. Figure 4 5 indicates that 
the distribution for each of the two lasers is a combi- 
nation of the distributions of each chaotic solitary laser, 
as shown in Figure 3a. The distribution consists of two 
maxima related to the refractoriness connected to pi and 
P2- For time intervals greater than both refractory times 
the distribution follows a Poisson distribution. Similar 
results are obtained for all combinations of p values ly- 
ing in the range of 1 to 1.5. Thus a statistical analysis 
reveals information about and differentiates between two 
sources operating with different parameters. 

Figure 4 5 indicates that although the two mutually 
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FIG. 5: Two lasers in a face-to-face configuration correspond- 
ing to a fully transparent mirror in the schematic of Figure 
1. The first laser operates with p = 1.1 and the second laser 
with p = lA and the coupling strength is a few % of the 
laser intensities. The probability for the time intervals be- 
tween spikes is presented for each one of the lasers. The two 
maxima correspond closely to the refractoriness of the solitary 
lasers operating with p = 1.1 and 1.4 (see Figure 3). 



interacting lasers with differing p values are not synchro- 
nized, the distribution of the intervals between spikes of 
each individual laser contains information about the pa- 
rameter p of the other laser. More precisely, from the 
measurement of the distribution of spiking time intervals 
containing the two refractory periods only and using its 
own parameter p, each laser can deduce the parameter p 
of the other laser. This mechanism may play an impor- 
tant role in neurobiology, where the statistical measure of 
the short intervals among spikes may reveal information 
about the individual state of the interacting neurons. 

For the third scenario of a partially transmitting 
mirror, each laser receives both self-coupling and mu- 
tual coupling signals, and the two lasers synchronize 
isochronally with zero time lag. The correlation coeffi- 
cient for the intensity traces averaged over 200 100 ns 
long time segments, exceeds 0.95 for p > 1.2 and is 
around 0.9 for small p, where low frequency fluctuations 
[23i] appear which are taken as partand are included in the 
calculation of the statistics. In the following we investi- 
gate whether the deviation from perfect synchronization 
is caused by a mismatch between the timing of the spikes, 
by the difference in their heights or by the background 
noise in our measurements. 

Figure 5a 6a shows the histogram for the mismatch 
between the timing of the spikes of the two lasers in the 
isochronal phase. The histogram indicates that the most 
probable time difference between the spikes of the two 
coupled lasers is zero since the average difference between 
the timing is less than 25 ps, our sampling rate. The 
width of the histogram, limited by our detection band- 
width, is extremely narrow (< 80 ps), and is not resolved 
by our detection system. 

We also examine the relative difference between the 
maximum intensities of two correlated spikes in the two 
lasers, and the histogram consisting of over 60, 000 pulses 
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FIG. 6: Two lasers in a zero-lag isochronal phase with a par- 
tially transmitting mirror in the schematic of Figure 1. (a) 
The histogram of the mismatch in the timing of spikes of the 
two lasers, and (b) the relative ratio between the maximum of 
the intensity of temporally correlated spikes of the two lasers. 



is presented in Figure 5b6b. Shown is the difference be- 
tween the maximum of each spike for each laser divided 
by the average of the spike amphtudes Aamp/Aveamp- 
This result indicates that the average relative difference 
between the maximum heights of temporally correlated 
spikes is ^ 10%. It is possible that the mismatch between 
the maximum intensities of these spikes is much lower, 
since the maximum intensity measured is very sensitive 
to the precise sampling time of the oscilloscope relative 
to the shape of a spike. From these observations we can- 
not definitively determine the source of the nonperfect 
synchronization, but it appears likely that it is not in the 
spike timing but rather in the spike amplitude. 

Phenomenon similar to our observations have been 
found in the communication of neurons, where imme- 
diately after the activation of an action potential it is 
more difficult to excite a second spike. Neural communi- 
cation has been documented to have many features (such 



as refractoriness, the repetitive form of the spiking pat- 
tern, synchronization between spatially separated neuron 
groups, the spike statistics [HjH, [2^) which are similar 
to those observed for the chaotic lasers. Our demon- 
stration that both the timing and maximum intensities 
of spikes are extremely well synchronized with zero time 
lag could have possible important implications for corre- 
sponding neural system [26, , ,27j . One of the fundamental 
problems in neuroscience, for example, is the question of 
how information is encoded in the neuronal spike trains 
[28| . Is the information contained in an individual spike 
form or in the interval between spikes, or is it the mean 
rate of spikes and timing which matter (29j? Tradition- 
ally it has been thought that most of the relevant in- 
formation was contained in the mean firing rate of the 
neuron. It is clear, however, that an approach based on 
a temporal average neglects all the information possibly 
contained in the exact timing of the spikes and the sta- 
tistical measure of the short intervals among spikes may 
reveal information about the individual state of the inter- 
acting neurons.. Recently more and more experimental 
evidence has accumulated which suggests that a straight- 
forward firing rate concept based on temporal averaging 
is too simplistic to describe neural information transfer. 
If, at each processing step, neurons had to wait and per- 
form a temporal average in order to read the message, 
the reaction time would be incompatibly long compared 
to experimental evidence. 



In conclusion, our results for chaotic lasers show that 
individual spiking laser units are able to generate irregu- 
lar spike patterns which become synchronized when two 
such units are coupled to each other, without any time de- 
lay, although the transmission time can be relatively long. 
Synchronization is maintained even on the time scale of 
individual spike widths. For chaotic lasers, transmission 
of information by the spiking pattern has been demon- 
strated, and the repetitive bar-code pattern we observe 
has features which are useful for communication appli- 
cations of these signals. The similarity of the lasers to 
neural systems is noted, and it is possible that complex 
neural systems can be effectively modeled by the much 
simpler and much more flexible and well-controlled ex- 
perimental environment of coupled laser systems. 



[1] H. Schuster, W. Just, Deterministic Chaos (Wiley VCH, 
2005) 

[2] P. Colet, R. Roy, Opt. Lett. 19, 2056 (1994) 

[3] G. D. VanWiggeren, R. Roy, Science 279, 1198 (1998). 

[4] A. Argyris, D. Syvridis, L. Larger, V. Annovazzi-Lodi, 

P. Colet, I. Fischer, J. Garcia-Ojalvo, C. R. Mirasso, L. 

Pesquera, K. A. Shore, Nature 438, 343 - 346 (2005). 
[5] E. Klein, R. Mislovaty, I. Kanter, W. Kinzel, Phys. Rev. 

E 72, 016214 (2005). 
[6] E. Klein, N. Gross, E. Kopelovitch, M. Rosenbluh, W. 

Kinzel, L. Khaykovich, I. Kanter, Phys. Rev. E 74, 



046201 (2006). 

[7] Jia-Ming Liu, How-Foo Chen, Shou Tang, IEEE J. Quan- 
tum Electronics 38, 1184 (2002). 

[8] S. Tang, J. M. Liu, IEEE Journal of Quantum Electronics 
39, 708 (2003). 

[9] T. Heil, I. Fischer, W. Elsasser, Phys. Rev. Lett. 86, 795 
(2001). 

[10] A. Hohl, A. Gavrielides, T. Erneux and V. Kovanis, Phys. 

Rev. A 59, 3941 (1999). 
[11] S. Sivaprakasam, E. M. Shahverdiev, P. S. Spencer, K. 

A. Shore, Phys. Rev. Lett. 87, 154101 (2001). 



5 



[12] A. Murakami J. Ohtsubo, Phys. Rev. A 65, 033826 

(2002). 

[13] R. Vicente, T. Perez, C.R. Mirasso, IEEE Journal of 

Quantum Electronic 38, 1197 (2002). 
[14] H. Erzgraber, D. Lenstra, B. Krauskopf, E. Wille, M. 

Peil, I. Fischer, W. Elsasser, Opt. Comm. 255, 286 

(2005). 

[15] E. Klein, N. Gross, M. Rosenbluh, W. Kinzel, L. 
Khaykovich, 1. Kanter, Phys. Rev. E 73, 066214 (2006). 

[16] I. Kanter, N. Gross, E. Klein, E. Kopelowitz, P. 
Yoskovits, L. Khaykovich, W. Kinzel, M. Rosenbluh, 
Phys. Rev. Lett. 98, 154101 (2007). 

[17] N. Gross, W. Kinzel, I. Kanter, M. Rosenbluh, L. 
Khaykovich, Optics Comm. 267, 464 (2006). 

[18] I. Fischer, R. Vicente, J.M. Buldu, M. Peil, C.R. Mi- 
rasso, M.C. Torrent, J. Garcia-Ojalvo, Phys. Rev. Lett. 
97, 123902 (2006). 

[19] A. Cho, Science 314, 37 (2006). 

[20] M. J. Berry, M. Meister, J. of Neuroscience 18, 2200 
(1998). 



[21] F. Ricke, D. Warland, R. de Ruyter van Steveninck, W. 

Bialek, Spikes - exploring the neural code, (MIT Press, 

Cambridge MA, 1996). 
[22] R. Lang, K. Kobayashi, IEEE J. Quantum Electron. QE- 

16, 347 (1980). 

[23] V. Alhers, U. Parlitz, W. Lauterborn, Phys. Rev. E 58, 
7208 (1998). 

[24] A. K. Engel, P. Koenig, A. K. Kreiter, W. Singer, Science 
252, 1177 (1991). 

[25] P. R. Roelfsema, A. K. Engel, P. Koenig, W. Singer, 
Nature 385, 157 (1997). 

[26] W. Gerstner, W. Kistler, Spiking neuron models (Cam- 
bridge University Press 2002). 

[27] M. Abeles, Firing rates and well-timed events, in Models 
of Neural Netwoks II, (Springer- Verlag, New- York 1994). 

[28] B. W. Connors, M. J. Gutnick, Trends Neuroscience 13, 
99 (1990). 

[29] P. R. Gray, Biophys J. 7, 759 (1967). 



